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ABSTRACT 

In a composite system of gravitationally coupled stellar and gaseous discs, we perform linear stability analysis 
for axisymmetric coplanar perturbations using the two-fluid formalism. The background stellar and gaseous discs 
are taken to be scale-free with all physical variables varying as powers of cylindrical radius r with compatible 
exponents. The unstable modes set in as neutral modes or stationary perturbation configurations with angular 
frequency uj = 0. The axisymmetrically stable range is bounded by two marginal stability curves derived for 
stationary perturbation configurations. By the gravitational coupling between the stellar and the gaseous disc 
components, one only needs to consider the parameter regime of the stellar disc. There exist two unstable regimes 
in general: the collapse regime corresponding to large-scale perturbations and the ring-fragmentation regime 
corresponding to short-wavelength perturbations. The composite system will collapse if it rotates too slowly 
and will succumb to ring-fragmentation instabilities if it rotates sufficiently fast. The overall stable range for 
axisymmetric perturbations is determined by a necessary D-criterion involving of the effective Mach number 
squared (i.e., the square ratio of the stellar disc rotation speed to the stellar velocity dispersion but scaled by 
a numerical factor). Different mass ratio S and sound speed ratio ij of the gaseous and stellar disc components 
will alter the overall stability. For spiral galaxies or circumnuclear discs, we further include the dynamical effect 
of a massive dark matter halo. As examples, astrophysical applications to disc galaxies, proto-stellar discs and 
circumnuclear discs are discussed. 

Subject headings: hydrodynamics — ISM: general — galaxies: kinematics and dynamics — galaxies: spiral — 
galaxies: structure — waves. 



1. INTRODUCTION 

Axisymmetric instabilities in models of disc galaxies have 
been investigated extensively in the last century (e.g., Safronov 
1960; Toomre 1964; Binney & Tremaine 1987; Berlin & Lin 
1996). For a single disc of either gaseous or stellar content, 
Safronov (1960) and Toomre (1964) originally introduced a di- 
mensionless Q parameter to determine the local stability condi- 
tion (i.e. Q> I) against axisymmetric ring-like disturbances in 
the usual Wentzel-Kramers-Brillouin-Jeffreys (WKB J) or tight- 
winding approximation. A more realistic model of a disc galaxy 
would involve both gas and stars as well as an unseen massive 
dark matter halo, all interacting among themselves through the 
mutual gravitation. ' Many theoretical investigations have been 
conducted along this track for a composite system of two cou- 
pled discs (Lin & Shu 1966; Kato 1972; Jog & Solomon 1984a, 
b; Berlin & Romeo 1988; Romeo 1992; Elmegreen 1995; Jog 
1996; Lou & Fan 1998b). In these earlier treatments, either 
an approach of combined distribution function and fluid or the 
formalism of two fluids have been adopted in a WKBJ modal 
analysis. While these model results were initially derived in 
various galactic contexts, they can be applied or adapted also 
to, with proper qualifications, relevant self-gravitating disc sys- 
tems including accretion discs in general, circumnuclear discs, 

' Magnetic field and cosmic-ray gas component are dynamically important on larg( 
therein) in the galactic gas disc of interstellar medium (ISM) but are not considere 



protostellar discs, planetary discs and so forth. 

The local WKBJ or tight-winding approximation has been 
proven to be a powerful technique in analysing disc wave dy- 
namics. Meanwhile, theorists have long been keenly inter- 
ested in a class of relatively simple disc models referred to as 
scale-free discs (Mestel 1963; Zang 1976; Lemos et al. 1991; 
Lynden-Bell & Lemos 1993; Syer & Tremaine 1996; Evans & 
Read 1998; Goodman & Evans 1999; Shu et al. 2000; Lou 
2002; Lou & Fan 2002; Lou & Shen 2003; Shen & Lou 2003; 
Lou & Zou 2004; Lou & Wu 2004). Scale-free discs, where aU 
pertinent physical variables (e.g., disc rotation speed, surface 
mass density, angular speed etc.) scale as powers of cylindri- 
cal radius r, have become one effective and simple vehicle to 
explore disc dynamics. Perhaps, the most familiar case is the 
so-called singular isothermal discs (SIDs) or Mestel discs with 
an isothermal equation of state and flat rotation curves (Mestel 
1963; Zang 1976; Goodman & Evans 1999; Shu et al. 2000; 
Lou 2002; Lou & Shen 2003; Lou & Zou 2004). In contrast to 
the usual WKBJ approximation for perturbations, perturbations 
in axisymmetric scale-free discs can be treated, in some cases, 
globally and exactly without the local restriction (i.e., valid only 
in the short-wavelength regime). It is therefore possible to de- 
rive global properties of perturbations. Using scale-free disc 

scales (Fan & Lou 1996; Lou & Fan 1998a, 2003; Lou & Zou 2004 and references 
. here for simplicity. 
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models, Lemos et al. (1991) and Syer & Tremaine (1996) both 
studied the axisymmetric stability problem for a single disc and 
found that instabilities first set in as neutral modes or stationary 
configurations with angular frequency uj = 0. 

The main motivation of this paper is to examine the global 
axisymmetric stability problem in a composite system of two 
gravitationally coupled scale-free discs. As a more general ex- 
tension to the previous two-SID analysis (Lou & Shen 2003; 
Shen & Lou 2003), we further consider a much broader class 
of rotation curves as well as the equation of state. This con- 
tribution gives an explicit proof that stationary configurations 
(w = 0) do mark the marginal stability in the two-fluid system, 
a cogent supplement to our recent investigation on stationary 
perturbation configurations (Shen & Lou 2004). 

2. TWO-FLUID FORMALISM 

As an expedient approximation, we treat both discs as razor- 
thin discs and use either superscripts or subscripts s and g to 
indicate physical variables with stellar and gaseous disc associ- 
ations, respectively. The large-scale coupling between the two 
discs is primarily caused by the mutual gravitational interac- 
tion. In the present formulation of large-scale perturbations, we 
ignore non-ideal diffusive effects such as viscosity, resistivity 
and thermal conduction, etc. It is then straightforward to write 
down the basic coplanar fluid equation set for the stellar disc in 
cylindrical coordinates (r, 9, z) within the z = plane, namely 



^ + ;^(^^V)+^-(EY): 







ot or 
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where S' is the surface mass density, is the radial bulk flow 
velocity, / is the specific angular momentum about the rota- 
tion axis along z-direction, 11* is the vertically integrated effec- 
tive (two-dimensional) pressure due to stellar velocity disper- 
sion and (p is the total gravitational potential. For the gaseous 
disc, we simply replace superscript or subscript s by g in the 
above three equations. The coupling of the two sets of fluid 
equations is due to the gravitational potential through the Pois- 
son integral 

-Gi:{r',il),tydr' 



(j)(r,9,t)= fd^p 



(4) 



[r'2 + r2-2rr'cos(?/'-6')]V2 ' 

where E = E* -i- is the total surface mass density of the com- 
posite disc system. 

The barotropic equation of state assumes the relation be- 
tween two-dimensional pressure and surface mass density in 
the form of 

U = Ki:", (5) 



:n/s:Eg-i, 



(6) 



where coefficients K >0 and n> are constant. This directly 
leads to the sound speed^ a defined by 

JEo 

which scales as oc Eg"' . The case of n = 1 corresponds to an 
isothermal sound speed a. 

From the basic fluid equations for stellar and gaseous discs 
above (the latter are not written out explicitly), we may derive 
the axisymmetric equilibrium background properties (Shen & 
Lou 2004). We presume that in the equilibrium background 
of axisymmetry, both rotation curves of the two discs scale as 
oc and both surface mass densities of the two discs scale 
as oc r~" where a and P are two constant exponents and the 
proportional coefficients are allowed to be different in general. 
By assuming the same power-law indices for the stellar and 
gaseous disc rotation curves (or equivalently the surface mass 
densities), it is possible to consistently construct a global ax- 
isymmetric background equilibrium for the composite disc sys- 
tem that meets the requirement of radial force balance at all 
radii [see equation (2) and the corresponding equation with the 
superscript s replaced by g] and satisfies the Poisson integral 
(4) simultaneously. 

The scale-free condition requires the following relationship 
among a, (3 and n (see Syer & Tremaine 1996 and Shen & Lou 
2004), namely 

1+4/3 



a = 1+2/3 and n- 



(7) 



1+2/? ■ 

Once the rotation curve is specified, aU other physical variables 
are simultaneously determined. 

With the knowledge of computing a gravitational potential 
arising from an axisymmetric power-law surface mass density 
for the background rotational equilibrium (Kalnajs 1971; Qian 
1992; Syer & Tremaine 1996), we can derive a self-consistent 
axisyimnetric background equilibrium surface mass densities 
as 



ys _ 
^0 - 



27rG(2/3Po)'-^"'^^(l + ^)' 



(8) 

where the coefficient Vo as a function of /3 involves F-functions 
r(-/3+l/2)F(/3) 



(9) 



2r(-/3+i)r(/?+i/2) 

and parameter (5 = Eq/Eq is the ratio of the surface mass den- 
sity of the gaseous disc to that of the stellar disc, A and D are 
two dimensionless parameters. We note that the value of 2f3P() 
falls within (0, oo) for the prescribed range^ of /3 G (-1/4, 1/2) 
and is equal to 1 when /? — > (i.e., the case of two gravitation- 
ally coupled SlDs). 

From the radial force balance of the background equilibrium, 
there also exists a relation such that 

(10) 



A?(0?+1): 



■AliDl+l). 



^ In a stellar disc, the velocity dispersion mimics the sound speed to some extent. 

' The valid range of /3 G (-1 /4, 1/2) is determined by (1) barotropic index n > for warm discs, (2) surface mass density exponent a < 2 so that the central point 
mass does not diverge, and (3) this /3 range is contained within a wider range of P & (-1/2, 1/2) (for a cold disc system) when the computed force arising from the 
background potential remains finite (Syer & Tremaine 1996). 
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where r/ = /A^ = /a^ is another handy dimensionless pa- 
rameter for sound speed ratio squared. We note that A is ac- 
tually the reduced'' effective sound speed [scaled by a factor 
(1 +2/3)'/^] and the parameter D = V/A where V is the re- 
duced disc rotation speed, is the effective Mach number [see 
eqns (11)-(12) later]. Condition (10) is very important in our 
analysis because the rotations of the two discs are not inde- 
pendent of each other but are dynamically coupled. It suffices 
to examine the parameter regime in either stellar or gaseous 
disc. In disc galaxies, the typical velocity dispersion in a stel- 
lar disc exceeds the sound speed in a gaseous disc, implying 
?7 > 1 so that inequaUty > holds. Therefore, the phys- 
ical requirement > absolutely guarantees > and it 
suffices to consider the stability problem in terms of Z)^ > to- 
gether with parameters 6 and r] for different values [note that 
= ?7(-Dj + 1)- 1]. With these explanations, we shall express 
other equilibrium physical variables in terms of parameters A 
and/). 

The specific angular momenta Jq and about the z-axis and 
the sound speeds a, and Ug in the two coupled equilibrium discs 
are expressed by 

fo=AsDy-^ , fo=AgDgr'-^, (11) 

a? = nK.mr' = A?/[(l +2(3^^] , 
al = nKgii:lr-'=Al/[(l+2pyf']. 

The disc angular rotation speed O = jo / and the epicyclic 
frequency k = [i2^/r)dir^fl)/dr]^^^ are similarly expressed in 
terms of A and D as 

=A,D,r-'-^ , Ks = [2(1 -/3)]i/'fi, , 
= A,Dgr-i-^ , Kg = [2(1 -/3)]i/2j7g , 
with d jo jdr = yk} j (20) to simplify later derivations. 

2.1. Linear Perturbation Equations 

For a composite system of two gravitationally coupled discs 
in rotational equilibrium with axisymmetry, we introduce small 
coplanar perturbations denoted by subscript 1 along relevant 
physical variables. The corresponding linearized perturbation 
equations can be derived from the basic nonUnear equations 
(l)-(4) as 



dYf, 1 d 



, S-^ dj\ 



^ + --(rI]^«S) + O.^ + -il^ = 



— + 0. — -2n,-=--^a,- + , 



dt 2a 



ul + ris 



dj\^_d_ 

80 



^0 



(14) 



for the stellar disc as well as their counterparts for the gaseous 
discs, together with the Poisson integral 



iir,0,t)= fd^ ^ 

J Jo \r'^ + r^-2rr' 



■Gm + J:l)r'dr' 



li/2 



(15) 



+ r^- 2rr' cos(V' - 0)] 

^ By the adjective "reduced", we refer to the part of a physical variable after 
V = Vr^^ , quantity V is referred to as the reduced disc rotation speed. 



relating the total gravitational potential perturbation (pi and the 
total surface mass density perturbation Si = Ej +Ej. 

Given a Fourier periodic component in the form of exp[i{u)t- 
mO)] for small perturbations in general, we write for coplanar 
perturbations in the stellar disc 

El = iJ,\r)exp[iiL0t-m9)] , 

u\ = U'ir)exp[iiu)t-me)] , (16) 

7i =J^ir)exp[iiujt-m9)] , 

as well as their counterparts in the gaseous disc, together with 
the total gravitational potential perturbation 

01 = V(r)e\p[i(u;t-m9)] , (17) 

where integer m is taken to be non-negative. For axisymmet- 
ric m = perturbations, we introduce Fourier decompositions 
in equations (14) -(15) for the stellar disc to deive 

iiOfj,' + -^iri:W') = , 
rdr 

2 

20, 

We do the same in parallel for the gaseous disc to derive 

r dr 

2 

20„ 



(18) 



(19) 



For the total gravitational potential perturbation, we simply 
have 

V(r)=6d^,' -G(M' + M^)cos(m^V^/ 



/() (r'2 + r2-2rr'cos'(/;)'/2 

Equations (18)-(20) are the basic coplanar perturbation equa- 
tions used for our axisymmetric stability analysis. 

2.2. Axisymmetric Stability Analysis 

For axisymmetric stability analysis with radial oscillations, 
we choose the Kalnajs potential-density pairs below because 
perturbations can be generally expanded in terms of such com- 
plete basis functions (Kalnajs 1971; Biimey & Tremaine 1987; 
Lemos et al. 1991; Lou & Shen 2003). Specifically, we take 

= cr*r"^/^exp(;^lnr) , fj,^ = r~^^^ exp(i^lnr) , 
y=-2^Gr(^-' + //)A/;„(0, 
where ^ is a 'wavenumber' characterizing the radial variation 
scale, C7* and cr* are two small real coefficients and the parame- 
ter function 

^ r(m/2 + /g/2 + 1 /4)r(m/2 - i^/2 + 1 /4) 

2r(m/2 + i^/2 + 3/4)r(m/2-iC/2 + 3/4) ^ ' 
removing the power-law radial dependence. For example, in the disc rotation speed 



(21) 
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is the Kalnajs function (Kalnajs 1971) that involves 
r-functions of complex arguments. Note that Nm is even 
in ^. It then suffices to consider only ^ > 0. 

Using the first mass conservations in equations (18) and (19) 
respectively, we infer U oc iujr^/'^'^^l^^^^ (see also Lou & Zou 
2004). By potential-density pair (21), equations (18) and (19) 
reduce to 



(uj^-Hi)U' = -G2U^ , 



(23) 



iL0^-H2)U«=-GiU' 
in the Umit of w — > 0, where parameter functions H\, H2, Gi 
and G2 are expUcitly defined by 

,(^'+l/4) 



■-27rG7VoSJ 



Gi = 2nGAfo^i^-^-^^^ > , 
r 

,(^2+1/4) 



(24) 



G2 =27rG7VoS«- 



>0 



The axisymmetric dispersion relation in the composite system 
follows from equation (23) 

uj'^-(Hi+H2)LJ^ + (HiH2-GiG2) = (25) 
in the limit of a; ^ 0, which is identical in form with earUer 
results obtained in the WKBJ regime (Jog & Solomon 1984a; 
Shen & Lou 2003). It would be of interest to note that the con- 
ditions for the stellar disc and the gaseous disc to be separately 
stable are Hi>Q and H2 > 0, respectively. It is also reminded 
here that in the familiar WKBJ regime, the dispersion relation 
in a single disc is 

uj^ = K^+ k^a^ - 27rGkSo , (26) 
where k is the radial wavenumber, whereas in the present global 
analysis as appUed to a single disc, we have in the limit of 
w->-0 



:.^,(i!lVl)«?_2.G^^^±i/l^Eo. 



(27) 

r 

If we replace M)(0 approximately by (^^+1/4)"'/^ in the 
asymptotic regime of ^ 1, we readily identify the correspon- 
dence between the effective wavenumber (^^+1 /4)'/^ and kr 
in the WKBJ limit of kr ^ 1 by directly comparing disper- 
sion relations (26) and (27). Physically, dispersion relation (27) 
is more generally apphcable beyond the WKBJ regime and is 
globally accurate only in the limit of w — > (Shu et al. 2000). 

Getting back to the composite system, there are two real 
roots^ of uj'^ of equation (25) 

ul = ^{(Hi+H2)±l(Hi+H2f-4(HiH2-GiG2)V^^} , (28) 

with uil root being always positive^ (Shen & Lou 2003). For 
the purpose of axisynnmetric stabiUty analysis, we only need to 
examine root, namely 

wl=^{iHi+H2)-[im+H2f-4iHiH2-GiG2)V^^} • (29) 

' One can show that the determinant of equation (25) A = (Hi — 7/2)^ +4Gi G2 
^lfHi+H2> 0, then u)J > 0; otherwise if Hi +H2 < 0, then at least one of Hi 



As the right-hand side of the above equation (29) is always real, 
axisymmetric instabihties set in as stationary perturbation con- 
figurations with ijl = and leads to the marginal stability con- 
dition 

HiH2 = GiG2 (30) 

that requires inequality H1+H2 > 0. This inequahty can be 
shown in a straightforward manner to be automatically satisfied 
if equation (30) holds true. Let us first write 



Hi=Fi-Gi and H2 = F2-G2, 
where Fi and F2 are explicitly defined by 



(^'+l/4) , 



a^>0, F2 = kI + 



(^'+1/4) , 



(31) 

>0. 
(32) 



It then follows from condition (30) of H1H2 = G1G2 that 

FiF2 = FiG2 + F2Gi . (33) 

As Fi, F2, Gi and G2 are aU positive, we immediately conclude 
that 

Fi -Gi > and F2-G2 > , (34) 
which finally leads to 

Hi +H2>0. (35) 

The composite disc system becomes inevitably unstable for 
H1+H2 < because of uj^ < 0. Else if Hi+H2>0 but with 
H\H2-G\G2 < 0, we again have < for instabilities. Only 
when Hi +H2 > and H1H2 - Gi G2 > at the same time can the 
composite disc system be stable against axisymmetric coplanar 
perturbations. This is an important necessary stability criterion 
for a composite system of two gravitationally coupled discs. In 
other words, once one disc is unstable by itself (i.e., < or 
7/2 < or both), the two-disc system must be unstable; even 
if the two discs are both separately stable, the composite disc 
system can still become unstable (i.e.. Hi >0 and H2>0 but 
//1//2-G1G2 <0). 

By inserting expressions of Hi, H2, Gi and G2 into the 
marginal stability condition (30) together with requirements 
of the background rotational equilibrium, we readily derive a 
quadratic equation in terms of y = D^, namely 



C2/ + Ciy + Co = 0, 
where the coefficients are explicitly defined by 
C2 =BonoV , 

(Ao + BoXHo-Bo) 



(36) 



Ci 



Co. 



(Bo-Ao)no + - 
(Ao + BoKHo+BoS) 



(l + S) 



- AqHo + 



+ iAo+Bof 



(A,+Bo)(no- 



(37) 



Bo) 



(1 + 5) 
iA)+Bo)ino+Bo6) 



(l + S) 

> is always true. 

and H2 is negative. It therefore follows that H1H2 — G1G2 < and hence i^J > 0. 



Axisymmetric stability 



5 



and 

A(o=e'+i/4, 

Bo(/?) = (l + 2/3)(-2 + 2/3), 
Ci(3) = (1 +2(3) /(2pVo), 

This quadratic equation (36) of y = can be readily proven 
to always have two real solutions (Shen & Lou 2004). Only 
the positive portion of Dj solutions can be regarded as phys- 
ically acceptable. Typically, there exist two different regimes 
bounded by the marginal stabihty curves that are unstable 
against axisymmetric perturbations, that is, the collapse regime 
for long-wavelength perturbations and the ring-fragmentation 
regime for short-wavelength perturbations. In contrast to the 
short-wavelength WKBJ approximation, the collapse regime is 
novel and exact. Systems with too fast a rotation parameter 
will fall into the ring-fragmentation regime (Safronov 1960; 
Toomre 1964; Syer & Tremaine 1996; Lou & Fan 1998a, b; 
Shu et al. 2000; Lou 2002; Lou & Shen 2003, 2004; Shen 
& Lou 2003), while those with too slow a parameter will 
fall into the collapse regime. Shown in Fig. 1 is an exam- 
ple of illustration with /3 = 1/4, ?y = 1 and an unconstrained 5. 
We note that cases with ?7 = 1 are essentially the same as those 
of a single disc (Shen & Lou 2004). The boundaries of the 
two regimes of instabilities shown in Fig. 1 vary with different 
parameters r] > I and S for chosen values of /3 that describes 
the entire scale-free radial profile of a composite disc system 
(i.e., disc rotation curves, surface mass densities and barotropic 
equation of state). Qualitatively, the increase of either r] and S 
will aggravate the ring-fragmentation instabiUty while suppress 
the large-scale collapse instabiUty (Shen & Lou 2003). While 
this can be directly seen from equation (36) in the relevant pa- 
rameter regimes, it can also be understood physically in terms 
of the dynamical couphng between the two discs through con- 
dition (10) constrained by scale-free conditions. For a larger rj, 
the reduced gas disc rotation speed Vg will exceed the reduced 
stellar disc rotation speed Vj by a larger margin and this tend to 
prevent an overall collapse of a composite disc system. In other 
words, a gaseous disc component with a relatively lower sound 
speed seems to prevent a collapse (Shen & Lou 2003). More 
details and examples can be found in Shen & Lou (2004). 

It is well known that a Q parameter can be defined to de- 
termine local axisymmetric stability of a single-disc system 
(Safronov 1960; Toomre 1964; Binney & Tremaine 1987). For 
a composite disc system, it has been attempted to introduce an 
effective Q^ff parameter (Elmegreen 1995; Jog 1996; Lou & 
Fan 1998b). As a result of straightforward numerical computa- 
tions, we have recently introduced a powerful D-criterion for 
axisymmetric stability of a composite system of two coupled 
SIDs (Shen & Lou 2003). It is natural to further generalize this 
D— criterion for axisymmetric stability of a composite system 
of two coupled barotropic discs by straightforward numerical 
computations as shown in Fig. 1. In the present case, it is more 



practical and simple to use parameter. 
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Fig. 1 . — One illustrating example for the two unstable regimes witli m = 0, 
/3 = 1/4 and rj = I. The collapse regime is at the lower-left corner, while the 
ring-fragmentation regime is at the upper-right corner In this special case of 
»7 = 1, parameter 5 can be arbitrary as can be seen from equation (37), namely, 
when r)= I the coefficients C2, Ci and Co turn out to be independent of <5. The 
vertical dash-dotted line is the location of where TIq = and C2 vanishes. 
The solid line and the dashed line bound the collapse regime and the ring frag- 
mentation regime, respectively. Only when falls within the range between 
the top of the collapse regime and the bottom of the ring fragmentation regime 
can a composite disc system become stable against all axisymmetric coplanar 
perturbations. 



2.3. Partial Disc Systems and Applications to Disc Galaxies 

From observations of more or less flat rotation curves of most 
disc galaxies, massive dark matter halos have been inferred to 
exist ubiquitously as long as the Newtonian gravity remains 
valid on galactic scales. If we naively attempt to relate theo- 
retical results obtained in Section 2.2 to a typical disc galaxy, 
we may take the simple isothermal equation of state as an ex- 
ample of illustration. The relevant parameters for a compos- 
ite SID system are then chosen as /3 = and = 50 km s~\ 
Vj = 220 km s"' , ^ = 0. 1 and rj = 50. Unfortunately, such a com- 
posite system of two coupled SIDs is inevitably unstable against 
ring-fragmentation perturbations because of a Z)j ~ 20, far ex- 
ceeding the maximum value necessary for being stable against 
ring fragmentations. This dilemma can be resolved by attribut- 
ing to an additional gravitational potential associated with an 
unseen massive dark matter halo. We refer to a composite disc 
system in association with an axisymmetric dark matter halo 
as a composite system of partial discs (e.g. Syer & Tremaine 
1996; Shu et al. 2000; Lou 2002; Lou & Fan 2002; Shen & Lou 
2003, 2004). In a simple treatment, the dynamical effect of a 
dark matter halo is modeled as only to contribute an axisym- 
metric gravitational potential $ in the background rotational 
equilibrium but not to respond to coplanar perturbations in the 
composite disc system. With ^7- = $ + ^, we conveniently in- 
troduce a dimensionless parameter = ip/ipr as the ratio of the 
potential arising from the composite disc system to that of the 
entire system including the presumed axisymmetric dark matter 
halo. The full-disc system corresponds toJ^=l (i.e. $ = 0) and 
the partial-disc system corresponds to < < 1 (i.e., $ j^O). 
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For a composite system of two gravitationally coupled partial 
discs, we follow the same procedure of analysing coplanar per- 
turbations in full discs to derive a similar quadratic equation of 
as the stationary dispersion relation. In a nutshell, we can 
simply replace all MoiO out theoretical results by J^MoiO 
to accomplish this generaUzation or extension. The introduc- 
tion of the ratio parameter T will significantly reduce both the 
ring-fragmentation regime and the collapse regime, as already 
can be seen from a comparative study of the WKBJ and global 
approaches (Shen & Lou 2003). 

For the purpose of illustrating the stabilizing effect of a 
partial-disc system, we simply take = 0.1 with other param- 
eters used earlier in this section. The minimum value of 
for unstable ring fragmentations now becomes ~ 650, far be- 
yond the actual value of Z), ~ 20 in a disc galaxy. Meanwhile, 
the collapse regime completely disappears. Therefore, a typical 
composite system of two coupled partial discd is fairly stable 
against axisynometric coplanar perturbations. 

3. DISCUSSION AND SUMMARY 

The main thrust of this investigation is to model linear copla- 
nar perturbations of axisymmetry (m = 0) in a composite system 
of two-fluid scale-free discs with one intended for a stellar disc 
and the other intended for a gaseous disc. The two discs are 
dynamically coupled through the mutual gravitational interac- 
tion. In order to include the dynamical effect of a massive dark 
matter halo with axisymmetry, we further describe a compos- 
ite system of two coupled partial discs (e.g. Syer & Tremaine 
1996; Shu et al. 2000; Lou 2002; Lou & Shen 2003; Lou & 
Zou 2004; Lou & Wu 2004). In a global perturbation analy- 
sis, we show that axisymmetric instabilities set in as station- 
ary perturbation configurations with a; = 0. The marginal 
stabiUty curves (characterized by the stationary configurations) 
deUneate two different unstable regimes, namely, the collapse 
regime for large-scale perturbations and the ring-fragmentation 
regime for short-wavelength perturbations. Apparently, the 
composite disc system becomes less stable than a single-disc 
system and can be unstable with the two discs being stable sep- 
arately (Lou & Fan 1998b). In our analysis, stationary perturba- 
tion configurations turn out to be more than just an alternative 
equilibrium state, especially in view of the stability properties. 

The basic results of this paper are generally apphcable to 
self-gravitating disc systems with or without axisymmetric dark 
matter halos. The two-fluid treatment contains more realistic 
elements than a single-disc formulation in the context of disc 
galaxies. In addition to astrophysical applications to disc galax- 
ies, the studies presented here can be valuable for exploring 
the dynamical evolution of protostellar discs and circunmuclear 
discs. 

In the context of a proto-stellar disc, it is the usual case to 
ignore the self-gravity effect. By considering the self-gravity 
of a composite disc system, our analysis indicates several qual- 
itative yet interesting results. For example, if the initial disc 
system rotates sufficiently fast, the ring fragmentation (see the 
upper-right part of Fig. 1) can occur at relatively small radial 



scales. By further non-axisymmetric fragmentations, these con- 
densed rings of materials may eventually become birthplaces 
of planets. On the other hand, if the initial disc system rotates 
sufficiently slow, then gravitational collapse can be induced by 
perturbations of relatively large radial scales (see the lower left 
comer of Fig. 1). Once such a perturbation develops in the 
background equilibrium disc, it grows rapidly and destabilizes 
the disc. Subsequently, the system undergoes global Jeans col- 
lapse to form a central young stellar object. Finally, if the initial 
disc system rotates in a regime stable against all axisymmetric 
perturbations (see Fig. 1), there might be two possibilities: (1) 
the composite disc system might become unstable caused by 
non-axisymmetric perturbations (not analyzed here) and (2) the 
disc rotation may be gradually slowed down by some braking 
mechanisms (e.g. magnetic field not included here and outflows 
or winds) and the disc eventually succumbs to a central collapse 
induced by large-scale perturbations. 

Likewise, in the context of a circumnuclear disc around the 
center of a galaxy, we can readily conceive similar physical 
processes in parallel. One important distinction is that a dark- 
matter halo should play an important dynamical role so that a 
formulation of a partial composite disc system would be more 
appropriate. Here, the ring fragmentation can be induced by rel- 
atively small-scale perturbations in a disc system of sufficiently 
fast rotation. Such a ring of relatively dense materials around 
the galactic center would be a natural birthplace for circumnu- 
clear starburst activities (e.g. Lou et al. 2001). Depending on 
the evolution history of a circumnuclear disc system, it may be 
stable initially and gradually lose angular momentum by gener- 
ating and damping spiral magnetohydrodynamic (MHD) den- 
sity waves (Lou et al. 2001). When the disc rotation becomes 
sufficiently slow, Jeans collapse induced by large-scale pertur- 
bations can set in to form a bulge or a supermassive black hole. 

In summary, our global analysis shows the possible presence 
of an evolution stage for a composite disc system against all 
axisymmetric coplanar perturbations. More importantly, we re- 
veal the parameter regime of ring fragmentation and the param- 
eter regime of large-scale collapse. Astrophysical apphcations 
are discussed in the contexts of disc galaxies, proto-stellar discs 
and circunmuclear disks. 
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